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1. Introduction 

RBC/UKQCD Collaboration has been performing a large scale simulation of the lattice QCD 
with dynamical domain-wall fermion (DWF) [|I|]. In this project, we can intensively access the 
CKM matrix, which includes b-quark physics. To treat the b-quarks, the Heavy Quark Effective 
Theory (HQET) is widely used. The lattice calculation with the HQET, however, has following 
difficulties (and solutions). 

1. The static propagator is too noisy. — This is basically because the static self-energy contains 
I /a power divergence. ALPHA collaboration investigated carefully this phenomena and 
introduced a modified static action which improves the signal to noise ratio [|3|]. The modi- 
fication can be achieved by replacing the link variable in the static action with the smeared 
one, which is obtained by the 3-step hyper-cubic blocking. Using this the power divergence 
contributions in the static self-energy are largely reduced. 

2. Non-perturbative matching with continuum is needed. — If we include 0{l/mh) correc- 
tion in the HQET formulation, the continuum limit cannot be reached by using perturbative 
matching factor because of power divergence [Q]. Possible ways of non-perturbative match- 
ing are Schrodinger functional scheme with step scaling technique and RI/MOM scheme. 

While the calculation can, in principle, be performed using the techniques described above, the 
actual implementation is not easy. There is an established way to apply the RI/MOM scheme for 
DW light quarks. But it has not been applied to the HQET successfully. 

As the first step of the project, the static approximation (lowest order of the HQET) is valuable 
and an important approach to the complete HQET. In the static limit, the perturbative matching 
procedure is justified. The perturbative matching at 0{as) in the static heavy and DW tight quark 
system was calculated without tink smearing in [||] and with link smearing in ||^ In this report, 
we present the calculation including the 0{a) improvement, whose effect cannot be neglected in 
the heavy quark system that we are considering here. 



2. Action setup 

We use the Iwasaki gluonic action and DW fermion with light quark mass niq for the light 
quark sector. The DWF has an optimized parameter M5 which is called "domain-wall height" and 
takes value < M5 < 2. In the calculation of the matching factor, it is assumed that the extension 
of the 5th dimension is infinity, which means the light quarks have exact chiral symmetry. For this 
sector, we do not carry out the link smearing. 

For heavy quark sector, we use the static approximation with link smearing: 



J^h{x,t) h{x,t) -W^ {x,t - \)h{x,t - I] 



(2.1) 



where h{x,t) is the effective heavy quark field and Wo(^,f) is the time-component of the smeared 
link variable. If Wq = Uo, which is the original gauge link, the action describes the one proposed 
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by Eichten and Hill [j8|]. We use the 3-step hyper-cubic blocked link for Wq with three parameters 
(«! , a2i 0:3). Possible parameter choices are 

(0.0,0.0,0.0) 

(ai,a2,«3 



(1.0,0.0,0.0) 
(0.75,0.6,0.3) 
(1.0,1.0,0.5) 



unsmeared [Wq = Uq) 
APE with a = 1 @ 
HYPl 
HYP2 B. 



(2.2) 



3. 0{a) in the static heavy and light quark system 

In this report we mainly focus on the 0{a) improvement of operators, and then we treat the 
matching factor between continuum HQET (CHQET) and lattice HQET (LHQET). Perturbative 
matching at one-loop between continuum QCD and CHQET was obtained by Eichten and Hill 
which we can use. 

Quark bilinear operator 



We consider the on-shell 0{a) improved static heavy (h) - light (q) quark bilinear 



O' 



CHQET 



Zr(l +bYtnqa) 



(1) 



(3.1) 



relating the CHQET operator O™^^^ on the left hand side and LHQET operators on the right hand 
side. Oy'' = hVq and O^' = hVy-Dq with T = {1, 7^, 75, 7^/5, CJ^v}- -Zr is the overall matching 
factor between CHQET and LHQET, cp and by are the 0{pa) and 0{mqa) improvement coeffi- 
cient, respectively. In this expression, we reduced the dimension 4 operators using the equation of 
motions of static heavy and massless light quarks 



Dq/j = 0, 



pq = 0. 



(3.2) 



The 0{pa) improvement of the heavy-light currents with clover Wilson light quarks was investi- 
gated using one-loop perturbation theory in non-relativistic QCD [O] and the static approximation 



[12]. They showed that the 0{pa) effects give a large correction to the B meson decay constant /g. 
In the light-light quark system, the existence of chiral symmetry guarantees the absence of 0{a) 
errors in the operators. For the case of the static heavy-light quark system, however, there are 0{a) 
effects even if we use chiral fermions for the light quarks. This was already found in the clover 
Wilson light quark with Wilson parameter r = (It is chirally symmetric, but there are doublers.) 



[12]. 



Now we consider the symmetries which the theory has. In addition to the chiral symmetry 
in the light quark sector, we have the heavy quark symmetry h —>■ e^'^i^i'''^'''h for the heavy quark 
sector. These symmetries guarantee that Zp is independent on T [|l^]. And also, cr = Gc, br = Gb 
with 7or7o = GV, in which c and b are independent on P. 

Four-quark operator 

2) which is relevant for the 



We consider the four-quark operator (AB 
(full) QCD operator is 



B mixing. Its 



(3.3) 
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where yj^ = J^Pl and also = J^Pr. At the one-loop level we need to take into account only the 
CHQET operator 

OL = 2[h^+Y^^0-¥^<l], (3.4) 

for matching between the CHQET and LHQET operators, where h^^\h^^^) is the particle (anti- 
particle) of the static quark. With the use of chiral fermions for the light quark, the on-shell 0{a) 
improved four-fermion operator can be written in 



^CHQET 



(3.5) 



where 



O^u = 2[/j(+)T^g][/j(-)T^(T4-A)'7]+4[/j(+)/ig][/j(-)p«(T^-A-)^] 

+2[U+^y^{riDi)q][h^-'>yf;q] +4[U+^PR{riDi)q0-^PLq], (3.6) 
On = 2p+^y^q0-^y^q]+4p+^PLq0-^PRq] 

+2[h'^+^Y^q]p-Y^q]+0^^PRqW-^PLq], (3.7) 

Zl is an overall matching factor, z|^^ is the 0{pa) improvement coefficient and z|"'^ is the 0{mqa) 
improvement coefficient. 

4. One-loop perturbative calculation of the 0(a) coefficients 

Quark bilinear operator 

We calculate the matching factor and the 0{a) improvement coefficients using one-loop per- 
turbation theory. The calculation is performed by comparing the light to heavy scattering amplitude 
between the CHEQT and LHQET. Now we consider the scattering amplitude with an initial light 
quark q carrying momentum p and a final heavy quark h carrying momentum k. In order to extract 
the on-shell 0{a) coefficients, the amplitude is expanded in the external quark momenta p and k 
around zero momentum and the light quark mass niq around zero mass. Since the momenta obey 



the equation of motions ( |3.2| ), / = and k^ = 0. In the puerturbative calculation, we choose the 
Feynman gauge and the UV divergences in the continuum calculation are regulated by dimensional 
regularization and we use the MS scheme for the renormalization. The IR divergences are regulated 
by inti^oducing the gluon mass A. 

The renormalized scattering amplitude for the CHQET at one-loop order can be written in the 
form 

m\OT\q{p)).ont = (l + ^Cf<L°2) {0?)o + ^CfS^}oUo?)o 

+ ^CF^}:-lm,a{0'^\, (4.1) 

where ( )o represents the tree level expectation value of the amplitude, Cf = {N^ — 1)/ (2A'^c) with 
number of color Nc, and 

^(o)__3, 5 M)__^ J,n) 4n 
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unsmeared 


APE 


HYPl 


HYP2 


5M 


12.979 


5.514 


4.910 


3.671 


e 


14.884 


1.429 


0.667 


-3.378 


eR = e- 5M 


1.906 


-4.085 


-4.243 


-7.049 



Table 1: Numerical values of dM, e and cr for each link smearing. 



In this expression n is the renormalization scale parameter. The scattering amplitude for the 
LHQET has the same IR divergence as in the continuum: 



{hik)\Jr\qip)) 



latt 



/ ,(0). 
latt j - 



wl)m,a{4%, (4.3) 



where wq - 
■^latt 



I-M5, 



3aX 



^ " ' -^latt 



4n 

3aX 



+ (4.4) 



The value of / was obtained in [[141]. Since we will use the fitting function ~ , = e — dM, 
the reduced value of e, is used [j8|]. The values are presented in Tab. [T[ d^^\ d^^^ and d^'"'^ are the 
finite parts of the vertex correction whose values are shown in Fig. [l[ After the matching we obtain 
the renormalized operator with 0{a) improvement 



where 



Zr 



l + -^Cf 

471 

An 



^ln(aV)+^ 



f + eR 



^CpGd^" 
An 



(4.5) 

(4.6) 
(4.7) 



Because of our use of DW light quarks we need the DW-specific factors ( 1 



0; 



'l-Cl-M^ 



7b are new results of this 



and Zh. = 1 + ^CfZw [jlj] in Eq. (^. The 0{a,a) coefficients Eq. 
calculation. 

Here we should briefly mention the \/a power divergence in the operator O^^ caused by the 
mixing with the lower dimensional operator O^^. Since this power divergence is already in the 
0{a) part, the total contribution is 0{a^) and we do not worry about it in taking the continuum 
limit. And also this 0{a^) effect contributes at 0{aj), which we can neglect in this one-loop 
calculation. This is quite different from the power divergence that appears in the 1 / mt, expansion: 
if the matching is done at Z-th loop, the power divergence of al^^/a remains. 

Four-quark operator 

The calculation for the four-quark operator can be done just by rearranging the quark bilinear 
results above. After the matching we obtain the 0{a) improved operator 



O 



CHQET 



1 - (wo)")"nZvv)"'ZL OL + zl^^aOND + Zr'{l-w^)mgaON 



(4.8) 
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Figure 1: Numerical values of d'^^\ d^^^ and (1 — m'q)c/^'"'' versus DW height M5 for each link smearing. 



As in the quark bilinear operator, we need the DW-specific factors in Eq. (4.8). The matching 
factor and 0{a) coefficients are 



7(0) 

7(1) 



1 + 



471 



/ 2 2^ 7 10 ,m c V 4 4 „ 



471 ' ^ 471 



(4.9) 
(4.10) 



where the constant v is the one-loop contribution from the diagram in which the gluon connects 
light and light, which was obtained in [0]. The constant c arises when the gluon connects two 
heavy lines and is given by c = cr. 

5. Discussion 



In Eqs. (|3JJ) and ( |3.5| ), we used the operators Op and Ond which contain covariant deriva- 
tives. These operators, however, can be written in the form: 

/jWr/.D^ = TGdo(h^^^rq) , 0ND = 2[h^+^y^q] (7o-7o) [h'^-^T^q], (5.1) 



where we have used the equations of motion ( |3.2| ). This form is quite convenient for taking the 
0{a) improvement in correlation functions. In the evaluation of the 2-point correlation function 

(a5)"^™P(0aJ)"^^(0)), where a[)"^ = h'^-'^YoYsq, we have 

(Af,-'™P(OA;,")^(0)) = {l+br{l-wl)m,a + CAaE^ind) {Ai-\t)A'^-^\o)) . (5.2) 

£^bind is the binding energy of static heavy and light quark, which is obtained in the correlator fitting. 
Therefore, in order to accomplish the 0{a) improvement, no further measurement is needed in the 
2-point correlation function. And also because the 0{mqa) part can be neglected due to its small 



size in many cases, we omit the 0{mqa) in the following discussion. Using the Eq. (5.2), we can 
evaluate the 0{a) improvement of B meson decay constant /g like 



/™P = (1 +CAa£'bind)/B- 



(5.3) 



For matrix element B parameter Bb and SU(3) breaking ratio t, we use the vacuum saturate 
approximation (VSA), and we obtain 



pimp 



VSA 



l+2(zi' 



CA)a£'bmdj Bb, 

1 +Z\^'a{E\,{nA(B,) — -£'bind(Brfj 



:-P ( 1 ^yfD.^p. £hind^«.)) ) ^ (5.4) 
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Now we roughly estimate these 0{asa) effect using the actual simulation data (j8 = 2.13, 
X r X L5 = 16^ X 32 X 16, M5 = 1.80, niuda = {0.01,0.02,0.03}, m,a = 0.0359) which ap- 
peared in [^. For this estimate the MF-improvement is taken into account. In this case, J^'' = 
3.48(APE), 6.41(HYP2) and £bmd ~ 0.6(APE), 0.5(HYP2). The coupling constant has the range 
a, 0. 15 — 0.35, conservatively. The conclusion is that the 0{asa) effect of /g is 3 — 8% (APE) , 5 — 
12%(HYP2), of is 9 -24% (APE), 15 - 36%(HYP2), and of is 3 - 8% (APE), 5 - 12%(HYP2). 
Using the assumption (£'bind(Bj) — ^bind(Brf)) ~ ("IBj — fnE^), the effect for E, is less than 2%. 

6. Summary 

We have presented a one-loop perturbative calculation of the 0{a) improvement coefficient 
for the static heavy - DW light quark system taking into account the hnk smearing in the heavy 
quark sector. Estimated 0{a) effect is not small in /g, and Bb, but is small in E, . While pertur- 
bative matching has large ambiguities and its own limitations, we deduce that this conclusions is 
not largely changed even in the non-perturbative matching. 

We thank all the member of the RBC and UKQCD Collaborations. 
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